Random resistor network model of minimal conductivity in graphene 
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Transport in undoped graphene is related to percolating current patterns in the networks of N- 
and P-type regions reflecting the strong bipolar charge density fluctuations. Transmissions of the 
P- N junctions, though small, are vital in establishing the macroscopic conductivity. We propose a 
random resistor network model to analyze scaling dependencies of the conductance on the doping 
and disorder, the quantum magnetoresistance and the corresponding dephasing rate. 
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Graphene - an atomic monolayer of graphite [l| - is a where 
gapless semiconductor with linear Q electron spectrum. 
The carrier density in a graphcne-bascd field effect tran- 
sistor (GraFET) can be varied continuously from P-type 
to iV-type. Experimental results [ij are surprising: the 
resistance per square never exceeds several fcfi's in con- 
trast to the pinch-off in conventional (gapfull) semicon- 
ductors. Such a behaviour of GraFET persists over a 
broad temperature range 10~ 2 K < T < 10 2 K. The 
conductivity reaches its minima when the gate-controlled 
carrier density, n g = 0, - "neutrality point" . Similar ob- 
servations were reported for the bilayer graphene 

Attempts to understand the finite value of the minimal 
conductivity in graphene have so far addressed the role of 
"chirality" of the Dirac-type quasiparticles and influence 
of various short-range defects on the quantum transport 
Q in a homogeneous graphene sheet. However, there 
is an emerging evidence [a H that charge density in 
GraFETs is "mesoscopically" inhomogeneous - probably, 
reflecting the fluctuations of the charge trapped in the un- 
derlying substrate, or on its surface. An inhomogeneity 
of charge density in a "charge-neutral" graphene sheet, 
n g = implies that it can be viewed as a checkerboard of 
N- and P-type doped regions separated by weakly con- 
ducting Q P-N junctions. In this Letter we propose a 
random resistor network model for such a system and use 
it to describe classical and quantum transport. 

Model of the random resistor network (RRN) is for- 
mulated on the square lattice with the lattice constant 
a and sites labelled by integers Sites 
with \i' — i\ > 1, \j' — j\ > 1 are not connected directly. 
Each pair of sites with \i' — i\ < 1, and \f — j\ < 1, is con- 
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nected by a link with the conductance 



G 



(i 



j) 



g [1 + /2; (la) 



+1 d) = 9 [1 



= G 



Here rjij 



is a random variable 



(-1) 



i+J 



19: 



m,j] /2; 

7 < 1. 



±!j {Vi,j)=P, {Vi,jVk,l) = SikSji, 



(lb) 

(lc) 

(Id) 




FIG. 1: Random resistor network (RRN) representation of a 
graphene sheet with nominally zero doping, n g = 0. 

The relation of model |T|) to graphene with inho- 
mogeneous charge density, n(r) = n g + 5n(r) is illus- 
trated on Fig. [T] Assume that the density fluctuations 
6n(r) are characterised by the length scale a so that 
(Sn(r)Sn(r + x)) = Sn 2 f(x/a); a may be determined by, 
e.g., the thickness of the insulating substrate. As long 
as a 2 Sn 3> 1, one can view the system as a combination 
of electron (N) and hole (P) puddles of the size of the 
order of a. Each puddle contains a large number of car- 
riers, which are characterised by the Fermi wavevector 
kp ~ \Jixb~n » a -1 . If i + j is even (odd), the site of 
the RRN corresponds to the A r (P)-puddle and is marked 
by red (blue) colour on Fig. [TJ The correspondence be- 
tween the puddles and the lattice sites encodes the fact 
that the observable conductivity is determined by the 
random links between the puddles rather than the local 
conductivity of a puddle. 

Each plaquette of the RRN, see the inset for Fig. [TJ 
has one and only one diagonal connection - either P-P 
or N-N link. This is described by Eqs. ((Taj) and flb| . If 
parameter p = 0, then P-P and N-N connections appear 
with equal probability, and p > (p < 0) describe the 
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electron (hole) doping. The boundary between the pud- 
dles has a finite, though small, transparency Q. This 
transparency is characterised by 7 <C 1 in Eq. (fTc|) and 
depicted by the square-lattice grid on Fig. Q] 

Scaling analysis - First, let us consider a RRN with 
p = 0, 7 = in Eq. ([T]), when two-colour random net- 
works at a bipartite lattice have peculiar geometrical fea- 
tures illustrated in Fig. [T] by a computer-generated sam- 
ple. Emerging patterns are typical for the percolation 
theory [9(. The RRN is critical, i.e. the geometry of 
the RRN is self-similar on all length scales: L x L net- 
work contains, typically, a larger cluster of one polarity 
(such as the red cluster in the middle part of Fig. [T]) 
which separates a pair of smaller blue clusters. Each of 
those clusters is in turn a shell for several smaller red 
clusters, etc. Such alternating cluster-embedding repre- 
sents a scale-invariant property of an infinite network. 
As a result, larger and larger parts of the network be- 
come excluded from the mono-colour percolation upon 
the increase of the size L. Therefore, the observable con- 
ductance G(L) decreases with the increase of L. The 
corresponding critical behaviour of its mean value is Q 



(G(L)) ~ (o/L)*g; x « 0.97. 



(2) 



This means that the conductivity is not defined. 

Finite p and 7 are relevant perturbations for the per- 
colation leading to a finite correlation length 7). At 
L ^S> £, the RRN is not critical and consists of indepen- 
dent patches of size £. Thus, the conductivity is finite: 



a(j,p) = (G(L ^ cx>)) ^ [a/^(p,j)] x g. 



(3) 



Because of the scale invariance, £ depends on p, 7 as 

£ (p, 7 0) ~ a\p\- v ; £ (p 0, 7) ~ a^. (4a) 

Its p-dependence is the dependence of the correlation 
length on deviations from the percolation threshold Q: 



4/3. 



(4b) 



To the best of our knowledge, the scaling of £ with 7 was 
not considered in literature, yet. We found 



/i = l/(h + x) « 0.37, 



(4c) 



where x is the coductance exponent, Eq. ([2]), and h = 
7/4 (see Ref. Q) is the exponent of the cluster outer 
perimeter V{L) ~ a(L/a) h for the cluster of the size L. 

To derive Eq. (|4c|) . consider the red cluster of the size 
L embedded into the conducting blue cluster. Small but 
finite 7 does not affect the criticality if the leak through 
the perimeter of this cluster V is much smaller than the 
conductance of the blue cluster, the latter is also of the 
order of G(L). One can use Eqs. ([2|) and (fTc|) to write 
this condition as (a/L) x > (L/a) hr y. This yields L < £ = 
0/7^, with [i given by Eq. (|4"c)) . For L > £, the leakage 



through the cluster boundary is efficient enough for the 
clusters to become distinguishable, i.e. RRN is uniform. 
The dependence of the conductivity on p, 7 7^ 0, is 



£ (p,7) ~ °7 '"/^(p/P*) > P*=7 



(5) 



where is a universal scaling function. The form of 
Eq. ([5]) is protected by the scale invariance. The numerics 
described below are well fit by the interpolation formula 



^(z) = (l + z 2 ) 
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(6) 



reproducing both z <C 1 and z 3> 1 asymptotic be- 
haviour. Substituting Eqs. ([5]), © into Eq. ^ we find 

o"(7,p) = ug^ a J^ ip/p*) , a = 2;^ « 0.36, (7) 

where u k 1.3 is a coefficient found from the fit of the 
numerical data. Equation ([7]) completely describes non- 
analytic dependence of the conductivity on both the leak- 
age parameter 7 <C 1 and adjustable doping p. Re- 
cently, an attempt was made in Ref. [HI to explain the 
graphene minimal conductivity using a mean field theory. 
We believe that Ref. ll| correctly describes the limit of 
high carrier density, where the fluctuations of local con- 
ductivity are small, but fails near the neutrality point 
where the percolation physics start to dominate. 

Numerical simulations of model ([T]) were performed on 
square lattices with L < 500a. Numerical evaluation of 
the network conductance has been done using the bus- 
bar boundary conditions at the two opposite edges along 
the direction of the current flow and the periodic bound- 
ary conditions in the orthogonal direction. Realizations 
of the RRN were generated in families spanning the in- 
terval p 6 [— 1; 1] of the network parameter, in steps of 
Ap = 0.0125. Each family was obtained from the reg- 
ular network at p = 1 (containing only blue links) by 
sequentially replacing L 2 /Ap blue links by red ones in 
randomly chosen plaquettes. Averaging of the data has 
been performed over 300 families. To minimise the finite- 
size effects we (i) identified 1 / L corrections by finite size 
scaling from L = 50a to L — 500a; and (ii) subtracted 
those corrections by means of numerical extrapolation. 
Numerical results summarised on Fig. [5] are in excellent 
agreement with the scaling formulas (O, ([7]). 

Sample-to-sample fluctuations - Geometry of the per- 
colation cluster is non-trivial. In particular, its lin- 
ear size depends on the particular realization and re- 
conecting only few bottleneck links may change it sub- 
stantially. Therefore, the conductance G(L) is a ran- 
dom quantity, which is characterised by its distribution 
function P(G) rather than only by its average. Scal- 
ing invariance of the critical cluster at L <C £(p, 7) 
constraints the functional form of this distribution to 
P{G) = (G)' 1 V{G/(G)), f¥{x)dx = 1, where P(x) is a 
universal function. It implies the variance 



(5G 2 (L«0) =u 2 {G{L^d)Y 



(8) 
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FIG. 2: (a) Conductivity <7 m in(7) = c(7,p = 0). Numerical 
data are represented by dots. The solid line is the best fit to 
the scaling law in Eq. (J7J). (b) Collapse of the conductivity 
data obtained for RNNs with various —1/2 < p < 1/2 and 
values of the parameter 7 onto a single curve, Eq. ©, rep- 
resented by solid line. For each value of 7 the corresponding 
values of the parameter p, were found using the by fitting 
a(p) to Eqs. @ and (J5]), see inset. 



where u-i ~ 1 is a universal numerical coefficient. For 
L>^, the central limit theorem is restored, and we find 

(SG 2 (L » 0) = <<5G 2 (6) {t/Lf = u 2 a 2 * (C/L) 2 . (9) 

Figure [3] (a) illustrates the fluctuations of RRN conduc- 
tance. One can see that not only value of the minimal 
conductance but also its position fluctuates from sample 
to sample. Figure [3] (b) shows an excellent agreement 
with the scaling form |[SJ for p = 0. 

The fluctuations © are self-averaging. Notice, how- 
ever, that because £ is large, those fluctuations can ex- 



ceed the universal conductance fluctuations 



7* lis 



which are not self-averaging as long as the phase coher- 
ence is preserved. The scale for the energy dependence of 
fluctuations © is determined by the properties of the in- 
dividual bottle-neck links in the cluster and significantly 
exceeds the Thouless energy for the UCF. This observa- 
tion is consistent with numerical findings of Ref. 13]. 

Quantum magnetoresistance - Up to now we ignored 
the quantum interference and interaction effects. Some, 
and most importantly those that give rise to the quan- 
tum magetoresistance 1J], can be analysed within the 
percolation cluster framework [l5| . The role of the weak 
localisation (WL) is to renormalize the conductance of 
each link in Eq. (JTJ) due to the interference of the ran- 
dom walk paths that go through the same link twice. 



-0.2 -0.1 0.1 0.2 



P 




100 



FIG. 3: (a) Conductance of several 100 x 100 samples with 7 = 
10 -4 as a function of p: The parameter p is gradually changed 
from - .2 to 0.2 by flipping a small number of randomly chosen 
links from P-P type to the N-N type. The fluctuations of 
the conductance are the strongest at p = 0. (b) R.m.s. of 
sample-to-sample fluctuations of RRN conductance for p = 
and various values of 7 as a function of system size, L. 



For WL magnetoresistance only long paths are impor- 
tant. Thus, one can assume that the return probability 
is independent of each particular links, so that it can be 
included into the renormalization of the conductance g. 

As usual , the evaluation of the WL correction con- 
sists in finding the classical return probability for the 
random walk, whereas the quantum mechanics deter- 
mines the relevant spatial scales and the pre-factors in 
this problem, which are the same as in the case of homo- 
geneous graphene [l6| ■ 

To find the probability for the returning path, one 
poses the problem of the diffusion on the lattice, Eq. |T]), 
where the link conductances G are replaced by the dif- 
fusances. Then, comparison with Eqs. ©, © yields the 
scale-dependent diffusion coefficient D(L) in terms of the 
observable diffusion constant at large distances, D^: 



D(L)~[Z/min(L,0} x D $ . 



(10) 



Our treatment of the WL differs from that of Ref. [16| 
only by the scale dependence of the diffusion constant. 
The magnetoconductance 5a(B) = cr(B) — cr(0) in not- 
so-low magnetic field B can be estimated as 
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Here t± is the inter- valley scattering time and 1/ry is the 
rate for intra-valley scatterings breaking certain symme- 
tries of the system, see Ref. [16( for more details. The 
phase relaxation time will be estimated below. 

If t& > i 2 /D^, Eq. (fTi) reduces to the results of 
Ref. [l6| ■ The specifics of the scale dependent diffusion 
(fT0|) are revealed in the opposite limit <C £ 2 /£)j: 
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(12) 



where T(z) is a universal function 17[ with the asymp- 



totic behaviour Y(z <C 1) ~ z 2 , and Y"(z ^> 1) « 
1 — z~ x / 2 . The overall magnetoresistance ([12]) is strongly 
suppressed in comparison with that for the homogeneous 
graphene, \5a\ <C e 2 /(n 2 h). 

The phase relaxation in ordinary disordered conduc- 
tors is dominated by the electron-electron interaction 
(l8| . The time defining the magnetoresistance curvature 
d 2 a/dB 2 \B^o is controlled by the dimensionless conduc- 
tance G{L) = 2irhG(L)/e 2 : 



Ufa ^ T/g(Lf); = (£>t ) 



1/2 



(13) 



where = £ is the dephasing length. As the all 
the geometric properties of the system in Eq. (|13|) are 
encoded into the scale dependences of G(L) and D(L), 
Eq. (fl~3"|) is valid even for diffusion along the critical clus- 
ter. We obtain from Eqs. (fT3j) . (|2j). (|3)). and (fTOj) 



L* = ^ (Tz/T) ; T £ = W^W? « £ 



-2(l+x) 



(14) 



where $(cc) is a scaling function with the asymptotic be- 
haviour $(z > 1) w i/i, $(z > 1) w z 1 /(2+2a;)_ Condi- 
tion T ~ determines the crossover from the critical to 
the normal diffusion: L^T > T c ) > f, L (T < T ? ) < f. 

Equations (H3J) and (Tl~2"|) predict interesting doping be- 
haviour of the magnetoresistance at fixed temperature. 
At p < 1, is small, T^> T and the usual logarithmic 
magnetoresistance occurs. For p — > 0, T£ < T, and the 
WL becomes suppressed as l/£ x . Such prediction seems 
to be consistent with the experiment of Ref. [l9[ . 

Relation to the experimental parameters - To relate 
Eq. (0 to the properties of graphene, we have to con- 
nect the parameters 7 and g with the physical param- 
eters of the sample. If the charge inhomogeneity is the 
dominant disorder in the graphene monolayer, we can 

estimate the conductance of N-N (P-P) connections be- 

2 

tween puddles as g ~ j-akp and conductance of the P-N 
junction separating puddles of the opposite polarity as 
ig ru zL( a k F ) l l 2 , see Eq. (2) of Ref. @], 7 ~ (ak F )~ 1/2 ~ 

(a 2 5ny 1/4 



According to Eq 

2 1 

CTmin = o-(-y,p = 0) ~ % (a 2 5n) 2 



(7) 



, 7 

we then estimate 



In conclusion, we constructed a random resistor net- 
work model, which adequately takes into account strong 
fluctuations of the local charge density and, thus, of the 
local conductivity of the mono- and bilayer graphene 
near neutrality point. This model describes inhomoge- 
neous current percolating through the system, giving rise 
to the scaling dependencies of the observable conductiv- 
ity on the doping and disorder. Quantum magnetoresis- 
tance and the sample-to-sample fluctuations are analysed 
within the model. 
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